Studies of Er ionization energy 

G. Gaigalas, Z. Rudzikas and T. Zalandauskas 

Institute of Theoretical Physics and Astronomy 
A. Gostauto 12, Vilnius 2600, Lithuania, 
t e-mail: tomas@itpa.lt 

February 2, 2008 



Abstract 

This work is aimed at the multiconfigurational Hartree-Fock calculations of the Er ionization 
energy. Authors have used the ATSP MCHF version in which there are new codes for calculation 
of spin- angular parts written on the basis of the methodology Gaigalas, Rudzikas and Froese 
Fischer 01 [5], based on the second quantization in coupled tensorial form, the angular momen- 
tum theory in 3 spaces (orbital, spin and quasispin) and graphical technique of spin-angular 
integrations. They allow the study of configurations with open /-shells without any restrictions 
and lead to fairly accurate values of spectroscopic data. 
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1 INTRODUCTION 



There is considerable interest in understanding the physics and chemistry of heavy atoms and 
ions. The main problem in investigation of such systems is their complexity, caused by the large 
number of electrons and the importance of relativistic effects. Therefore the detailed description 
of heavy atoms and ions requires the correct treatment of correlation effects as well as relativistic 
description of the system. The correlation effects in the theory of many-electron atoms are treated 
mainly by the two methods: configuration superposition (CI) and multiconfigurational Hartree- 
Fock (MCHF). Relativistic effects are usually included adding the relativistic corrections to the 
non-relativistic Hamiltonian (for example, the Breit-Pauli approximation), or using relativistic two 
component wave functions and relativistic Hamiltonian (for example, the Breit approximation) pQ. 
The complete and rigorous treatment of both relativistic and correlation effects for heavy atoms 
and ions is, unfortunately, practically outside of today's computational possibilities. 

Nevertheless, we expect that in case of Er relativistic effects as well as correlation effects between 
the electrons of 'inner shells' (core-core correlations) are the same for the neutral atom and ion. 
We assume that these effects (corresponding energies) cancel each other in calculation of ionization 
energy (IE). Only correlation between 'outer' (valence) electrons gives significant contribution to 
IE. And then it may be possible to get quite accurate values of ionization energies taking into 
account of them by MCHF approach. 

This work is aimed at checking this assumption. For this purpose we perform MCHF calculations 
using the ATSP MCHF J2J |3] version in which there are new codes for calculation of spin-angular 
parts written on the basis of the methodology Gaigalas, Rudzikas and Froese Fischer |U|!)|, based 
on the second quantization in coupled tonsorial form, the angular momentum theory in 3 spaces 
(orbital, spin and quasispin) and graphical technique of spin-angular integrations. They allow the 
study of configurations with open /-shells without any restrictions. 

2 APPROACH 

Ionization energy we define as IE = Ei on — E g , where E g and Ei on are the ground state energies 
of neutral and singly ionized atom correspondingly. The ground state of neutral Er is 

ls 2 2s 2 2p 6 3s 2 3p 6 3(i 10 4s 2 4p 6 4(i 10 5s 2 5/4/ 12 6s 2 3 H = [Xe]4f 12 6s 2 3 H (1) 

and that of singly ionized Er 

ls 2 2s 2 2p 6 3s 2 3p e 3d 10 4s 2 4p 6 4d w 5s 2 5p 6 4f 12 6s 1 4 H = [Xe]4/ 12 6s 1 4 H. (2) 

2.1 Hamiltonian 

We calculate the ground state energies making use of the Hamiltonian 
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H = H-NonRel + T~Lr 



elCor ; 



(3) 



where 7i NonRel is the usual non-relativistic Hamiltonian and TLReiCor stands for relativistic correc- 
tions. In our calculations we separate relativistic corrections into following three parts: 



7~(-RelCor = 7~(-Sh + 'Hmp + 7~LoO 

with mass-polarization correction given by the Hamiltonian 
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orbit-orbit term given by 
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and the remaining part of relativistic corrections 

Ti-sh = Hmc + Tim + Hdi + Ti-ssc, 
consisting of the mass correction term 

Kmc = -^r£- 



2 AT 

Pi, 
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(8) 



i=l 



as well as the contact interactions, described by the one- and two-body Darwin terms TCdi and 
H.D2- They are 
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i=l 



and 7Yd2 = — vra 2 E $ 

«7=1 



(9) 



The last addant in Eq. ((7J) stands for the spin-spin contact term 



7~Lssc 



8ira 
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E (si-Sj-)5(rij-) 



(10) 
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The corrections presented in Eq.(|%]) are (except mass polarization of the order of square of 
fine structure constant. They enable us to make a study of contribution of the main relativistic 
corrections to the calculations of ionization energy. Let us also mention that the Hamiltonian is 
presented in atomic units. 
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2.2 MCHF 



For calculation of ionization energy we used MCHF method. In this approach, the wave function is 
expressed as a linear combination of configuration state functions (CSFs) which are antisymmetrized 
products of one-electron spin-orbitals. A set of orbitals, or active set (AS), determines the set of 
all possible CSFs or the complete active space (CAS) for MCHF calculation. The size of the latter 
grows rapidly with the number of electrons and also with the size of the orbital AS. Most MCHF 
expansions are therefore limited to a restricted active space (RAS) [2|. No 'relaxation' effects were 
included. 

3 RAS construction 

Large scale systematic MCHF calculations of IE of lanthanides has not been done yet. Therefore, 
following the methodology of it is important to investigate the structure of Er and Er + ground 
configurations, to impose the core and valence shells and to evaluate valence-valence (VV), core- 
valence (CV) and core-core (CC) correlations. 

It is always a question when we can assume that a subshell is a part of the core, and when it should 
be treated as a valence shell. The answer is not trivial even for boron like ions, and in our case it 
is even more complicated because of complexity of configurations of Er and Er + , and our attempt 
is to take care of the correlation effects that do not cancel each other between ion and atom. 
Because we treat IE in non-relativistic approach, and in the neighbourhood of the ground level 
there are no levels with the same L, S, J values, methodics based on the consideration of energy 
spectra described in [2] could not be straightforward adapted to impose core and valence shells in 
our case. 

Therefore in this chapter we will study some possibilities of RAS construction. 

3.1 HF calculations 

First insight into the structure of Er and Er + ground states we can get from the Hartree-Fock (HF) 
calculations. The resultant ground state energies and mean distances of nl radial functions are 
presented in TABLE I. 

Resultant energies are in agreement with those presented in The important note is that 6s 
function is much more remote frome the nucleus than the ones of 5s, 5p and 4/. And the open 4/ 
shell is closer to the nucleus than the 5s and bp. Therefore, we have a difficulty in treatment of 
'outer' electrons: usually as outer (valence) shells the open ones are considered, but sometimes the 
closed shells are included too [2]. For light atoms these shells are spartially 'outer'. 

3.2 CORE I 

In this case we use the core 



4 



TABLE I. Results of HF calculations. Values of mean distance from the 
nucleus <r> and energies of ground states (in a. u.) 



Function 


<r> jtv 




Is 


0.022 


0.022 


2s 


0.094 


0.094 


2p 


0.080 


0.080 


3s 


0.242 


0.242 


3p 


0.232 


0.232 


3d 


0.205 


0.205 


4s 


0.545 


0.545 


4p 


0.557 


0.557 


M 


0.588 


0.588 


5s 


1.371 


1.385 


bp 


1.563 


1.565 


4/ 


0.754 


0.754 


6 s 


4.630 


4.093 


Energy: 


-12498.1528 


-12497.9809 



/ = [Xe] X S 

and the 4/, 6s we treat as valence shells. The 4/ shell we treat as valence shell because it is 
open and 6s because the corresponding radial function in much more remote from the nucleus than 
others. This approach is close to the advices given in [2]. 

The basis for the MCHF expansion was formed using the CSF's of configurations made of single 
and double (S, D) excitations from the valence shells to some destination set. There were two types 
of destination sets used: 

a = {5d, 5/, 5ss 6p, 6d} , (11) 

b = a + {6/, 6g, 6h, 7s, 7p, 7d} . (12) 

Further on we denote the basis as core with subscript of destination set. For example, I a denotes 
the basis, consisting of CSF's of configurations, made by S, D excitations from 4/ 12 6s 2 for Er and 
4/ 12 6s 1 for Er + to the destination set 'a' and cores [Xe]. The numbers of CSFs in the bases are 
presented in TABLE II. 

The weight for the main CSF was found to be 0.977 for I a (and similar for If,). This value is close 
to that (0.949) found by CI method [7]. The mean distances of radial functions from the nucleus 
are found to be smaller than for HF calculations. For example <r>4j = 0.752 a.u. for I a (0.748 
a.u. for lb) and <r>6 S = 4.550 a.u. for I a (4.534 a.u. for if,). 



5 



3.3 CORES II, III 



In this case, we treat as valence shell only 6s, because of its spatial location. We expect this 
strategy to be more efficient for the calculations of 6s ionization energy because as we can see from 
HF calculations mean distance of 4/ radial functions is not much different for Er and Er + . As a 
cores we use 

II. [Xe]4f 12 with not fixed term 
and 

III. [Xe]4/ 12 with fixed term 3 H. 

There were five types of destination sets used with these cores: 
(jllj) and ()12j) as for core I and three more 

c = b + {7f, 7g, 7h, 7i, 8s, 8p, Sd} , (13) 
d = c + {8f, 8g, 8h, 8i, 8k, 9s, 9p, 9d} , (14) 

e = d + {9/, 9g, 9h, 9i, 9k, 91, 10s, lOp, Wd} . (15) 

As we can see from TABLE II, the basis formed with the same destination sets is the biggest for 
the core I, the medium for core II and the smallest for core III. 

The weights of main CSFs in MCHF expansions are about 0.960 - 0.980 for all bases with cores 
II, III. The mean distance from the nucleus for 6s radial function is greater than the one obtained 
from HF calculations but smaller than obtained using bases with core I. For example, <r>g s = 
4.560 a.u. for III a , 4.564 a.u. for Ii7 M , e . 

4 6s IONIZATION ENERGY 

The results of MCHF calculations are presented in TABLE II. 

For the bases obtained using II, III cores in non-relativistic approach, we get that increasing 
destination set the value of IE decreases until certain value (in our case 5.792 eV). This value 
should be considered as the result of 'frozen core' method. The result that for certain core using 
smaller (for example 'a') destination set and correspondingly smaller basis we obtain the IE value 
closer to experimental one is treated by us as casual, because in smaller destination set (basis) we 
account for smaller part of valence correlations or (and) take into account it with different precision 
for Er and Er + . 

Adding the relativistic corrections raise the value of IE up to 5.877 eV. This value is fairly close to 
the experimental one 6.108 eV (the error is less than 4%). 
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TABLE II. Results of MCHF calculations. Numbers of CSFs (NCSF) and 

values of IE (in eV.) 



Basis 


NCSF (Er) 


NCSF (Er+) 


lENonRel 


lERel 


la 


2838 


1 1 by 


o.ooo 


o. i oy 


h 


12811 


12054 


5.572 


5.807 


Ha 


236 


8 


5.895 


6.640 


He 


2600 


23 


5.793 


5.877 


Ih 


5565 


32 


5.793 


5.874 


He 


10347 


43 


5.792 


5.877 


IHa 


70 


4 


5.896 


6.073 


Ilh 


272 


7 


5.796 


5.856 


Hie 


733 


11 


5.792 


5.876 


aid 


1569 


15 


5.792 


5.873 


Hie 


2938 


20 


5.792 


5.877 


Non relativistic HF 




4.677 


7 


CI 






5.077 


7. 


CIq with Davidson Q correction [S] 


5.250 


7. 


Estimated 


5.58 




7. 


Experiment 


6.108 




m 
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Using core I the corresponding basis is 10 times larger than the one of II (or about 40 times 
larger than formed using core III). Nevertheless, the results obtained using these bases are much 
worse. For example, in non-relativistic approach IE value for the If, basis is 5.572 eV, whereas the 
corresponding value for the If, is 5.796 eV. It is due to the fact that basis formed using destination 
set 'b' for the core I is not enough to account for the correlation effects of 4/-electrons, which, 
represented in full, cancels between Er and Er + . Relativistic corrections improve the value of IE, 
but it is still worse than those obtained using cores II, III. Of course, the results obtained using 
core I could be improved using larger bases (and should overtake the ones obtained using II, III 
cores), but at present it is unreachable because of our computational resources. 
So, we recomend to use core III for the calculations of IE. In this case it is possible to get quite 
good value of IE taking into account VV (6s-6s) correlations only. The inclusion of CV4/ and CC4/ 
correlations requires much more computational resources and with the present ones doesn't lead to 
the better result. 

And finally, let us compare our results with the previous calculations of IE, where correlation effects 
were treated using CI method 0. 

The authors used well-temped Gaussian type functions (GTF's) and augmented diffuse p-, d- and 
/-functions to describe the 6s-6s, 6s-4/ and 6s-5d electron correlations as well as 8g and 7h to 
describe the angular correlation effects. The total number of GTF's was (29s, 25p, 22d, 17 /, 8g, 7h). 
Since authors could not perform full single- and single-double excitation configuration interaction 
(SDCI), important CSF's were selected by performing the second order perturbation calculations, 
reducing number of CSF's. They performed so-called natural orbital iterations to obtain compact 
CI expansions for the ground state as well as for the ionized state. The error of the correlation 
energies due to unselected CSF's was estimated to be 0.001 a.u. and the error in the IE was 
estimated to be less than 0.05 eV. As one can see from the TABLE II, the result of CI calculation 
was 5.077 eV, and the one with Davidson correction 8 CIq 5.250 eV. The result with estimated 
relativistic effects is 5.58 eV. 

Comparing the results [7j with ours, (5.792 eV for non-relativistic value of IE), we can see that it 
is possible to obtain much better value of IE only accounting for the valence (6s-6s) correlation, 
whereas incomplete inclusion of correlation effects of inner shells just disimproves results. 

5 CONCLUSION 

For the calculations of Er ionization energy the most efficient strategy is to use the MCHF ex- 
pansions with frozen core [Ae]4/ 12 3 if and single, double excitations from (6s). CV4J and CC4J 
correlations are not important in this case. 

Our results on erbium are more accurate than data found using the Davidson Ci+Q method [S] 
and adding the relativistic corrections [7j . 
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The results obtained show that if the correlation effects of inner shells cancel each other between 
atom and ion, then it is possible to get quite accurate data by MCHF method accounting for 
correlation effects of outer electrons only. And this assumption takes place in the case of Er 
ionization energy. 

The results demonstrate the ability of the approach by Gaigalas et aLjlJEj to obtain fairly accurate 
data for heavy atoms and ions, having open /-shells. 
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